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Dynamics of Non Classically Reproducible Entanglement 
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We investigate when the quantum correlations of a bipartite system, under the influence of en- 
vironments with memory, are not reproducible with certainty by a classical local hidden variable 
model. To this purpose, we compare the dynamics of a Bell inequality with that of entanglement, 
as measured by concurrence. We flnd time regions when Bell inequality is not violated even in cor- 
respondence to high values of concurrence (up to ~ 0.8). We also suggest that these results may be 
observed by adopting a modification of a recent experimental optical setup. These findings indicate 
that even highly entangled systems cannot be exploited with certainty in contexts where the non 
classical reproducibility of quantum correlations is required. 
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In contexts such as quantum computation and quan- 
tum information, working protocols require for their im- 
plementation systems that present the peculiar quantum 
correlations characterized by entanglement P, Q ■ How- 
ever, unavoidable interaction of realistic systems with 
their environment gives rise to an increase of mixedness 
of the state of the systems and to a decrease of the de- 
gree of entanglement with time. Entanglement can even 
disappear completely at finite time (entanglement sudden 
death) . This motivates the interest in considering con- 
ditions and methods that maintain the systems entangled 
as long as possible. Among the conditions that effectively 
increase the entanglement usefulness time there are the 
use of non-Markovian environments, where entanglement 
revivals are possible 0, [1] , of quantum Zeno effect ,6] , or 
of structured environments, which can give rise to entan- 
glement trapping 0, Q . 

It has been however shown that there exist entangled 
bipartite mixed states whose correlations can be repro- 
duced by a local hidden variable model, that is by clas- 
sical systems Q, alt houg h they may still display "some 
hidden non-locality" [lO|. This indicates that, for mixed 
states, which are in practice the ones always encountered, 
a given value of entanglement by itself does not imply 
that their correlations cannot be classically reproduced 
with certainty. Hereafter, we refer to quantum correla- 
tions that are certainly non reproducible by a classical 
local model as inherently nonlocal correlations (INLCs). 
Because of the necessity of INLCs for device- independent 
and security-proof quantum key distribution protocols 
[m . [l^ and their relevance for quantum computation, 
it appears crucial to have indicators of their presence. 
One of such indicators is obviously given, for bipartite 
systems, by a Bell function and the presence of INLCs 
unambiguously identified by it violating a Bell inequality 
[isl 14 1 . From another side, the question is to deter- 



mine when quantum traits, linked to entanglement and 
suitable for quantum computation, are not classically re- 
producible with certainty. 



The aim of this paper is thus to consider one of the 
cases where the time when the entanglement is present 
can be extended and to compare it with the time regions 
when INLCs are present. The system we shall consider 
is that of two qubits in environments with memory (non- 
Markovian). These systems present in general entangle- 
ment revivals and one expects that, with an appropriate 
choice of parameters, also revivals of INLCs occur. Fi- 
nally, we shall examine if the conditions when this hap- 
pens can be obtained within the current experimental 
technologies. 

We take two independent non-causally connected 
qubits, each interacting with a distinct, but identical 
bosonic reservoir at zero temperature. Let {|0),|1)} be 
the basis states of the qubit and the operator O a pseudo- 
spin observable with eigenvalues ±1, defined as O = O-cr, 
where O = (sin ^ cos 0, sin sin 0, cos 0) is the unit vec- 
tor indicating a direction in the pseudo-spin space and 
(T — (ci, (72, fa) the Pauli matrices vector. The ex- 
pression of O in terms of the basis states is O{0, (p) = 
cos6l(|l)(l| - |0)(0|) + sin6l(e*'^|l)(0| + e-*"*|0)(l|). The 
Clauser-Horne-Shimony-Holt (CHSH) form of the Bell 
function associated to the two-qubit state p for the oper- 
ator O is Q 

B{p) = |(Oi(0i, 0i)O2(^2, h)) - (Ol(01, 0l)O2(02, 02))| 
+ (Oi(0;, 0i)O2(02, <^2)) + (Ol(0'i, 0l)O2(ei, 02)), 

(1) 

where (0102) = Ti{pOi02} is the correlation func- 
tion, with the index S = 1,2 referring to the 5'-th 
qubit. If, given the state p, a set of angles {(j)i,(j)2} and 
{0i,6[,92,0'2} exists such that the CHSH-Bell inequality 
B{p) < 2 is violated, the correlations cannot be simu- 
lated by any classical local model and are nonlocal. Such 
a set of angles always exists for pure entangled states but 
generally not for mixed states [l5| . 

The Bell function at time t is obtained from the two- 
qubit state p{t). In our system this can be determined, 
for any initial state, by the knowledge of the single-qubit 
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dynamics, whose exact solution is available when the 
reservoir is at zero temperature and has a memory (non- 
Markovian) In this case, the single qubit-reservoir 
evolution is represented by the quantum map, known as 
amplitude decay channel 



|0s)® |0i?) -> |05)®|0fl). 



|0i 



® |0k> + ^^|0s) ® (2) 



where jOg), jig) and |0_r), are respectively the states 
of the qubit S and of the reservoir i?, and p is the decay 
probability given by [lij 
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with d = ^/2^X^^}? . In Eq. ([3]), A defines the spectral 
width of the coupling and is connected to the reservoir 
correlation time t^j by A ~ r^^^, F represents the qubit 
excited state decay rate in the Markovian limit of flat 
spectrum and is linked to the system (qubit) relaxation 
time T5 by r « ^. We consider here the strong cou- 
pling regime defined by A/F < 2, where the reservoir 
correlation time is larger than the qubit relaxation time 
and memory effects become relevant. 

The two-qubit dynamics can be analyzed in a sim- 
ple way for any initial two-qubit state [^]; however, here 
we shall restrict the analysis to the pure Bell-like initial 
states /5$(0) = |$)('I>| and p*(0) = |*)(^'|, where 

1$) = a|01) +e**/3|10), |*) = a|00) + e**/3|lf), (4) 

with a,/3 non-negative reals and 0^+0^ = 1. These 
states belong to the class of X states whose density ma- 
trix px^ in the standard basis B = {|1) = |11),|2) = 
|10),|3) = |01),|4) = |00)}, has an X structure. This 
structure has been shown to persist during the time evo- 
lution due to the map of Eq. ([2]) and has the form [3] 



Px{t) 
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P22{t) 
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In our system the initially pure states ^$(0) = |<1')(<1'|, 
p^{0) — become, during the time evolution, the 

mixed states p<i,(t), p^(t) whose matrix elements as ex- 
plicit functions of time t and ofprobability amplitude a 
have been previously reported [J]. 

Since our goal is to find the existence of regions where 
B > 2, it is strategic to maximize, by an appropriate 
choice of the angles, the Bell function B{p) of Eq. ([T]). For 
any state px{t) the standard procedure [15] determines 
these angles as (f)2} = {{k + k')TT — S14 — S23, [k — 
k')Tr — Si4 -\- S23}, where k,k' are integer numbers and 
<5i4, S23 respectively the initial phases of pi4{t), p23{t), 
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FIG. 1: Maximum of Bell function, Bmax, in terms of the 
initial degree of entanglement (represented by a^) and of the 
dimensionless time Ft for the initial Bell-like states |$) = 
a|01) + P\W) (upper plot) and 1*) = a|00) + /3|11) (lower 
plot) at A/r = 10""^. The parts of the graphs emerging after 
the first collapse under the classical threshold value 2, "Bell 
islands", correspond to a return of INLCs. 
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in addition {6*1, 6'i, 6(2(0, 6*^(0} = {0,f ,arctan 
6*2 (t)}. With this choice the maximum of B is 

B^,^{px{t)) = 2^/p2(^) + Q2(i)^ (6) 

where 



r{t) = Pll{t) + P4i{t) ~ p22{t) 
Q{t) = 2(|pi4(t)| + |p23(i)|). 



P33{t), 



(7) 



The expression for B^i^^[px{t)) given by Eq. ^ coin- 
cides with the one which would be obtained using the 
formal Horodecki criterion (l7| . 

In Fig. [T] the regions where -Bmax > 2 are reported as 
function of time and of the initial degree of entanglement 
(represented by a^). The figure clearly displays regions 
of revivals of CHSH-Bell inequality violations (Bell Is- 
lands). These regions correspond to a return, after finite 
intervals during which Bmax < 2, of INLCs at space- 
like distances via local qubit-environment interaction. To 
quantify entanglement, we adopt the concurrence C [3], 
being C — Q for non-entangled states and C = 1 for max- 
imally entangled states. The time behavior of i?max and 
of C is plotted, for the same values of the parameters. 
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FIG. 2: Maximum of Bell function, Bmax — 2 (solid line) and 
concurrence C (dashed line) as a function of the dimension- 
less time Ft for initial Bell-like states |$) = (|01) -I- y2|10))%/3 
(upper plot) and 1*) = (|00) + V2\ll))/V3 (lower plot). The 
value of is fixed at 1/3 so that dark intervals of entan- 
glement (C = 0) occur for the state |*I'). The coupling pa- 
rameter is fixed so that X/T — 10"'* representing strong non- 
Markovian effects. 



in Fig. m The plot evidences the appearance of regions 
of entanglement (C > 0) where however the CHSH-Bell 
inequality is not violated (Bmax < 2). In these regions it 
is therefore not possible to say with certainty that quan- 
tum correlations are not reproducible by a classical local 
model From the figures one sees that Bmax reaches 
for the first time the classical threshold value Bmax = 2 
rather earlier than the first entanglement disappearance 
(C = 0). This behavior confirms what already found 
previously in comparing the dynamics of Bell inequality 
violation with entanglement decay for two qubits sub- 
jected to local decoherence in the Markovian limit at zero 
[l9| | and finite temperature [20|]. Successively, Smax re- 
mains below 2 until C reaches again a threshold value. 
In particular, for the initial state \^), we find that this 
threshold value, for any time, is given by 2al3/{l + a'^(3'^); 
the maximum value of this concurrence threshold is 0.8 
and it occurs when the state is initially maximally entan- 
gled (/3 = 1/ \/2) , while in the case represented in Fig. [2] 
its value is « 0.77. For the initial state \'^), the ana- 
lytical expression of the corresponding threshold value 
of concurrence is more complex and we do not report 
it; its maximum value is ~ 0.60 and it is obtained for 
f3 « 0.79, while for the case considered in Fig. [2] it takes 
the value « 0.59. We have then times of high values of 
concurrence when the CHSH-Bell inequality is not vio- 
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FIG. 3: Comparison between the maximum of Bell func- 
tion Bmax and the fixed angles function Ba^ in terms of the 
quantum channel decay probability p, for initial Bell states 
1$) = (|01)±110))/V2 (upper plot) and|*) = (|00)±|11))/^ 
(lower plot). 

lated. We also point out that the peaks of -Bmax decay 
much more quickly than the ones of C. 

In the analysis above, the revivals of INLCs have been 
analyzed in terms of Bmaxj which is obtained by using the 
time dependent angles 92{t) and 92{t). We however will 
show that, by suitably fixing these angles, the relative 
difference between the Bell function at fixed angles, Ba^, 
and the maximum Bmax, is in practice negligible when 
both are above the classical threshold B = 2. To this 
aim, we fix the angles so that Bfi^ coincides with B^i^^ 
a,t t — 0. For both the initial states /5$(0) = |$)($| and 
/5*(0) = |*)(^'| this happens for 62(0) = tt - 61^(0) = 
arctan(2Q;/3), which leads to 

BnMt)) = -j=J==[l + 4a' 13'- {2 + Aa'f3')p\, 

Bfi.(p*(t)) = ^=J==[l + 4a'f3'-Af3'{a' + l)p 

+4(3'p']. (8) 

These two quantities, for maximally entangled initial 
states, are plotted in Fig. [3] with the corresponding 
Smax(p*(0) and -Bmax(/5*(i)) as functions of the decay 
probability p. By inspection, for both states, Bgx ~ ^max 
above the classical threshold B — 2, the relative error 
Sr = (Bmax " -Bflx)/-Bmax falfiug in the range < i5r < 
10^'^. A realization of the dynamical Bell test that exper- 
imentally reveals "Bell islands" is thus possible, for any 
practical purpose, just by following Bfix, with the prac- 
tical advantage of not changing the angles with time. 

Conditions for observing the presence of entanglement 
revivals are within the range of current experimental ca- 
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pabilities as encountered in solid state devices or in cav- 
ity quantum electrodynamics [i, 21, S^l- Conditions re- 
quired for observing revivals of INLCs are more extreme 
and presumably, for the moment, non obtainable within 
these contexts. However, a recent optical experiment, ex- 
ploiting a Sagnac-like interferometer that simulates the 
amplitude decay channel with the decay probabili ty p 
linked to the rotation angle of an half- wave plate [23j . 
has made possible to follow the two-qubit entanglement 
dynamics for a form of p corresponding to Markovian 
reservoirs. As a matter of fact, this experimental proce- 
dure is valid for any p, independently from the values of 
parameters, so that by taking the explicit form of Eq. ([3]) 
one can in principle obtain, with the same setup, the time 
evolution of concurrence also for non-Markovian environ- 
ments, that is also for the dashed curves plotted in Figl^] 
In this optical context, the qubit states |0), |1) are coded 
by the H (horizontal) and V (vertical) polarizations of a 
photon, while the environment states are represented by 
two different momentum modes of the photon and pure 
entangled photon pairs are repeatedly generated for each 
value of the decay probability p [2^ . 

Here we suggest that, with appropriate modifications, 
the same experimental setup may also be used to follow 
the revivals of INLCs, overcoming the problem of reach- 
ing, in the previous contexts, the required physical condi- 
tions. This realization of the dynamical Bell test could be 
accomplished by substituting, in the optical setup previ- 
ously described, the final detectors with a standard Bell 
analyzer The orientation angle of each polarizer 

of the Bell analyzer, 6'p, corresponds to the half angle 
9 of the unit vector O, 9p = 6/2 [ij. The photon po- 
larization measurements have to be performed with the 
appropriate settings of the polarizer angles. In particu- 
lar, the initial states of the qubits with a = (5 = l/\/2 
are coded into |$+) = {\HV) + \VH))/^/2 or \^+) = 
i\HH) + \ VV))/V2 and the polarizers must correspond- 
ingly be set at the standard angles {^pi, ^pi, ^p2, ^p2} = 
{0,45°, 22.5°, 67.5°}[l3; while the initial states with 
a ^ -(3 ^ 1/V2 are |$_) = {\HV) - \VH))/y/2 or 
|*_) = {\HH) - \VV))lV2 and the standard angles are 
{6pi,d'^i,0p2,e[,^) = {0, 135°, 22.5°, 67.5°}. The realiza- 
tion of this Bell test would check the p-dependence of 
Bfix given by Eq. ([8]) and plotted in Fig. [3l yielding the 
observation of the "Bell Islands" with p given by Eq. ([3]). 

In conclusion, we have studied the dynamics of CHSH- 
Bell inequality violations for two independent qubits, 
each embedded in an environment with memory, and 
compared it with the dynamics of entanglement, as mea- 
sured by concurrence. We have individuated the time 
regions when the quantum correlations of the two qubits 
are not reproducible with certainty by a classical local 
hidden variable model. We have shown that there exist 
time regions when the CHSH-Bell inequality is not vio- 
lated even in correspondence to high values of entangle- 
ment (up to C ~ 0.8). We have also suggested that both 



dynamical behaviors here determined could be observed 
by adopting a modification of a recent optical experi- 
mental setup. The results here found indicate that the 
entanglement, even for rather high values of concurrence, 
cannot be exploited with certainty in those circumstances 
where the non classical reproducibility of quantum cor- 
relations is essential, as, e.g., in quantum cryptography. 
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